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INTRODUCTION 


GPS Overview 

The NAVSTAR GPS is a space-based, navigation satellite system that will provide 
the user with precise three-dimensional position and velocity information as well as 
coordinated GPS time continuously, regardless of location, altitude, weather or time of 
day [1]. GPS is separated into three major segments: 1) the space segment, which is 
comprised of the earth-orbiting satellite constellation, 2) the control segment, which 
monitors the health and orbits of the satellites and 3) the user segment, which is 
comprised of all of the air, land, sea and spaced-based users with GPS receivers. 

Space Segment 

The space segment, the constellation of earth-orbiting satellites, will consist of 21 
satellites plus three active spares, when fully operational in 1994. They will be 
distributed uniformly in six subsyncronous circular orbital planes providing a minimum 
of four visible satellites at any one time anywhere on earth. These orbital planes are 
inclined 55° with respect to the equatorial plane and separated longitudinally by 60°. The 
satellites orbit at an altitude of 10898 nm and with a period of 12 hours. 

Signals from the satellites are transmitted continuously on two L-band frequencies 
designated L, and U, where L x is 1575.42 MHz and Uis 1227.60 MHz. The L, 
frequency is modulated with both the coarse/acquisition (C/A) code and precision (P) 
code, while the I ^ frequency is modulated only by P code. Hie C/A code has a frequency 
of 1.023 MHz and a wavelength of 300 m, with a period of one ms. The P code has a 
frequency of 10.23 MHz and a wavelength of 30 m, with a period of one week. The 
structure and composition of the P code signal is classified by the U.S. Department of 
Defense. In addition, both the L, and Lj frequencies are further modulated by a 
navigation message which contains GPS time, GPS almanac data, satellite ephemeris 
data, satellite health, atmospheric propagation correction data as well as any other 
information needed by the GPS receiver. 

Control Segment 

The control segment monitors, and corrects if required, the satellite broadcast signals 
to ensure a pre-defined accuracy level. In addition, the control segment is responsible for 
monitoring and controlling the orbits of the satellites, maintaining the GPS time and 
uplinking necessary information to the satellites three times a day. The control segment, 
when fully operational, will consist of a Master Control Station (MCS), located at Falcon 
AFB, Colorado Springs, CO, five monitor stations, one located at the MCS and the other 
four on the following islands; Hawaii, Kwajalein, Diego Garcia and Ascension and three 
uplink antennas. The monitor stations allow simultaneous tracking of the complete 
constellation and relay orbital information, GPS time and any other necessary 
information to the MCS. The ranging information, acquired by the monitoring stations, is 
then processed for use in satellite orbit determination and systematic error elimination. 
The MCS then calculates corrections that are uplinked to the satellites via the three uplink 
antennas. v 

User Segment 

The user segment is intended for both military and civilian users of the system. In 
order to use GPS, a receiver is required. The GPS receiver consists of an antenna to 
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capture the GPS signals, an amplifier to boost the power level of the received signals and 
a digital computer to process the information contained within the signals. 

The GPS receiver selects and measures a minimum of four independent satellite 
signals in order to calculate a three-dimensional position fix. In addition, three- 
dimensional velocity and GPS time is calculated. The GPS design specifications require 
the calculated position to be accurate to within 15 m Spherical Error Probable (SEP), the 
calculated velocity to be accurate to within 0.1 m/s (la) and the GPS time to be within 
100 ns (la) [2], Since the GPS information is output in the World Geodetic System, 
1984 (WGS-84) Earth-Centered, Earth-Fixed (ECEF) coordinate system, military and 
civilian position data can be standardized on a worldwide basis. 

Position Solution 

The GPS receiver position fix is accomplished by means of passive tri-lateration. 
Since the positions of the satellites are known at all times, the GPS receiver position is 
determined by measuring the distance between the receiver and a minimum of four 
satellites. The satellite atomic cesium clocks are kept synchronized to the GPS time by 
the MCS. The less accurate clocks contained within the GPS receivers are not. 
Therefore, the measurement of the signal’s time of travel will be in error by an amount 
equal to the difference between the atomic time standard maintained by the satellite and 
the time maintained by the GPS receiver. This error produces an inaccurate range 
measurement between the GPS receiver and each satellite it is receiving, known as 
“pseudorange”. In addition to the above time difference error, a satellite clock error as 
well as ionospheric and tropospheric delays are present. However, the effect of the 
satellite clock error is negligible for the typical navigation solution. In addition, this error 
is indistinguishable from the ionospheric and tropospheric delays. Actual offsets of the 
satellite clocks are approximated by polynomials in time and transmitted as part of the 
satellite signal to the user for the correction of the measured pseudorange. The 
ionospheric and tropospheric delays can be calculated on the basis of ionospheric and 
tropospheric models. Therefore, for each pseudorange measurement, an equation can be 
derived that relates the measurement to the satellite position and the unknown quantities 
of GPS receiver position and GPS receiver clock error. Since four unknown quantities 
exist, a minimum of four pseudorange measurements must be taken and solved 
simultaneously to determine the three-dimensional GPS receiver position. 
Mathematically, the pseudorange is defined as [3]: 

p, = V(X - * )’ + (Y - y , ) 2 + (Z - z, ) 2 - c(dT) (1) 


where: 

Pi is the pseudorange measurement to satellite i 

X,Y,Z is the GPS receiver position in the ECEF coordinate system 

Xi,y ir Zi is the position of satellite i in the ECEF coordinate system 

c is the speed of light (299,792,458 m/s) 

dT is the GPS receiver clock error 

The effect of the GPS receiver clock error and the way in which its measurement is used 
to calculate the true position of the GPS receiver is shown in Figure 1. 
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GPS Receiver Clock Error And Calculation Of True Position 


This figure illustrates that the pseudorange radii do not converge at a single point but 
instead enclose the shaded triangular area. However, if the range value equivalent to 
c(dT) is added to (or subtracted from, if the case warrants) the pseudorange, then the radii 
converge to a single point, the GPS receiver position. Figure 1 represents only a two- 
dimensional case in which all three of the satellites and the GPS receiver lie in the same 
plane. However, the same logic applies to the three-dimensional case with four satellites. 

Geometric Effects 

GPS receiver clock error is not the only contributor to the position error. The relative 
geometry of the satellites being tracked and the GPS receiver also affects the position 
accuracy. Therefore, in order to determine the accuracy available from the satellites 
being tracked as a function of their geometry, the Dilution of Precision (DOP) values 
must be calculated [1]. 

The DOP is composed of Position Dilution of Precision (PDOP), which reflects the 
dilution of precision in three-dimensional position; Horizontal Dilution of Precision 
(HDOP), which reflects the dilution of precision in the two horizontal dimensions; 
Vertical Dilution of Precision (VDOP), which reflects the dilution of precision in the 
vertical dimension; and Time Dilution of Precision (TDOP), which reflects the dilution of 
precision in time, i.e., in the estimate of the range error due to the GPS receiver clock 
error. Geometric Dilution of Precision (GDOP) is a composite value composed of both 
the position and time DOP. Since the magnitude of the range error due to the GPS 
receiver clock error is multiplied by the DOP values to obtain the overall position 
accuracy, small DOP values are desirable [1]. The most frequently used measure of 
geometric performance is PDOP, which when multiplied by the range error due to the 
GPS receiver clock error determines the position error. Figure 2 relates satellite geometry 
to PDOP values. 
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GOOD (LOW) PDOP 


POOR (HIGH) PDOP 


Figure 2 

PDOP As A Function Of Satellite Geometry 


DGPS Concept 

DGPS represents a straight-forward method to significantly improve the accuracy of 
GPS. The general principle of DGPS is that by having a GPS receiver at a pre-surveyed 
location, the true values of the location are compared against the measured values, and 
the resulting pseudorange and pseudorange-rate difference corrections are sent real-time, 
to be applied to the airborne GPS receiver. This real-time uplink of the pseudorange and 
pseudorange-rate difference corrections enables the airborne GPS receiver to calculate its 
position relatively free of measurement errors due to ionospheric delay, tropospheric 
delay, ephemeris uncertainties, satellite clock error and selective availability (SA), 
intentional degradation of the C/A signal. SA is accomplished by the degradation of the 
broadcast ephemeris data (GPS orbits and clock offsets) and by the systematic 
destabilization of the GPS oscillator. The latter generates an irregularly changing error in 
the pseudorange and pseudorange-rate values. All of these errors are either eliminated or 
significantly reduced since they are common to both GPS receivers. 

DGPS research at the NASA Ames Research Center was initiated in the early 1980’s 
when a C/A code, single channel sequencing DGPS system was developed and flight 
tested on a NASA SH-3G helicopter [4-7]. The objective of the tests were to evaluate the 
use of DGPS to support helicopter terminal approach operations. Final approach 
positioning accuracy was 5.2 m (mean) ±8.0 m (2a) laterally and -7.7 m (mean) ±7.0 m 
(2a) vertically. The use of barometric and radar altimeters to enhance the vertical axis 
positioning accuracy yielded 5.0 m (mean) ±6.0 m (2a) vertically and 5.0 m (mean) ±4.0 
m (2a) vertically, respectively [5,7]. The time lags due to the satellite sequencing were 
too great for real-time guidance and navigation. However, the positioning accuracy 
attained during these tests demonstrated the potential of DGPS for terminal approach 
operations. 

Recent fixed-wing DGPS terminal approach and landing flight tests conducted at the 
NASA Ames Research Center yielded positioning accuracy of 0.1 m (mean) ±1.8 m (2a) 
laterally and -0.8 m (mean) ±6.6 m (2a) vertically, using P code GPS receivers [8]. Thus, 
a DGPS based guidance and navigation system has the potential to significantly enhance 
the unique capabilities of the helicopter, as well as improving pilot situational awareness. 
Specific helicopter missions that would benefit from such a system are: 

1. Low altitude operations 

2. Remote area landing 
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3. Search and rescue operations 

4. Inter-city operations 

The above missions could be carried out regardless of weather conditions, time of day or 
location 

The use of C/A code DGPS in real-time to provide high accuracy, precision 
navigation and guidance for helicopter approaches to landing is the subject of this 
research. The objective is to evaluate, through flight test, C/A code DGPS positioning 
accuracy during helicopter approaches to landing. 


5 


TEST EQUIPMENT AND FACILITIES 


Equipment Description 

The operational layout of the facilities for the flight tests is shown in Figure 3, which 
shows the various components involved. These include the NASA UH-60 helicopter 
with the Ashtech GPS receiver operating onboard, the ground-based Ashtech GPS 
receiver and antenna located at a pre-surveyed test location and the laser tracker used to 
provide the true position and velocity of the aircraft during the flight tests. 



AIRBORNE GPS 
RECEIVER 



Figure 3 

Flight Test Facilities 


Airborne System 

The test aircraft is a modified UH-60 Black Hawk helicopter located at the NASA 
Ames Research Center (See Figure 4). During the flight tests, two pilots are required to 
fly the aircraft and operate the aircraft systems. In addition, a research engineer is 
required to operate the airborne test equipment and to coordinate operations between the 
airborne and ground-based systems. The airborne test equipment consists of a Ashtech 
Model XII 12 channel C/A code GPS receiver, a Maxon/Ashtech SM 3010 VHF 
telemetry uplink receiver, a Litton LN-93 ring laser gyro Inertial Navigation System 
(INS) (used only for post-processing data collection), a 486 processor based Data 
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Acquisition/Navigation Computer (DAC) and an interface to the aircraft's standard 
approach guidance instruments (See Figure 5). In addition, a Tandy 102 laptop computer 
is used to control both the data acquisition and navigation functions. The GPS antenna is 
located at the rear of the aircraft on top of the vertical stabilizer, adjacent to the tail rotor. 
The telemetry uplink antenna is located on the bottom of the aircraft approximately 
midway between the main and tail landing gear. A laser reflector is mounted on both the 
right and left sponsons for position tracking during the approaches to landing. 

The DAC collects GPS position data at a rate of 2 Hz via a RS-232 connection and 
INS Euler angle data at a rate of 64 Hz via a 1553B data bus. In addition, a precise GPS 
time pulse is sent to the DAC at a rate of 1 Hz via a digital line. GPS time differs from 
Universal Coordinated Time (UTC), by not adjusting for leap seconds at periodic end-of- 
year intervals. Such adjustments in GPS time would disrupt the continuous availability 
of the satellites for navigation purposes. This time pulse is used to time-tag the INS data 
to that of the GPS time, via the IRIG-B time, which is derived from GPS time and 
broadcast from the test site. The telemetry uplink receiver is connected to the GPS 
receiver by a RS-232 connection and receives the pseudorange and pseudorange-rate 
difference corrections at a rate of approximately 0.5 Hz. The computer-derived lateral 
and vertical deviations are output to the approach guidance instruments at a rate of 2 Hz 
via a RS-232 connection. All airborne data collected from each approach is stored on a 
removable 44 MB Syquest hard drive for post-flight analysis. 


GPS 



Figure 4 

NASA UH-60 Helicopter 


GPS 



OPERATOR’S 

TERMINAL 


Figure 5 

Airborne System 
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The DGPS based guidance is displayed in the cockpit (using standard Instrument 
Landing System (ILS) type localizer and glideslope instruments) using a special interface 
unit that converts computer-derived lateral and vertical deviations from a nominal 
localizer and glideslope to simulated ILS localizer and glideslope signals, respectively 
[9]. 

Ground-Based Re ference System 

The ground-based test equipment consists of an identical Ashtech Model XII GPS 
receiver, a Maxon/ Ashtech SM 3010 VHF telemetry uplink transmitter and a laser tracker 
(See Figure 6). A research engineer monitors the ground-based test equipment and is in 
two-way radio contact with the aircraft at all times. The telemetry uplink transmitter is 
connected to the GPS receiver by a RS-232 connection and broadcasts the pseudorange 
and pseudorange-rate difference corrections at a rate of approximately 0.5 Hz. The GPS 
receiver antenna is permanently positioned at a pre-surveyed test location that allows for 
unobstructive viewing of the satellites. The telemetry uplink antenna is positioned so as 
to have line-of-sight with the aircraft throughout the flight test pattern. The laser tracking 
data is time-tagged with the IRIG-B time prior to being recorded. 


GPS 

ANTENNA 


UPLINK 

ANTENNA 



Figure 6 

Ground-Based Reference System 


Laser Tracking System 

All flight tests were conducted at Crows Landing NAS, Crows Landing, CA, located 
approximately 50 miles east of Moffett Field NAS. The NASA test facility at Crows 
Landing NAS includes two Nike X-band monopulse radar trackers and one precision 
NiYag laser tracker. The laser tracker provides precise aircraft range, azimuth and 
elevation and is used to provide the reference or truth data by tracking the reflectors 
mounted on the aircraft. Laser range accuracy is nominally ±0.3 m out to approximately 
9 km; azimuth and elevation accuracy are nominally ±0.2 mrad. During the terminal part 
of the approach this corresponds to position errors of less than 0.5 m in each axis. Note 
that these accuracies are of the same order of magnitude as the expected DGPS position 
accuracies. The laser tracker was calibrated each morning prior to a test flight. Refer to 
Appendix A for a geographic description of the test facilities at Crows Landing NAS. 
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DGPS BASED GUIDANCE SOLUTION 


The airborne GPS receiver calculates its position (m) and velocity (m/s) in the ECEF 
coordinate system at a rate of 2 Hz. In this reference frame, the origin is located at the 
earth’s center of mass, the X E axis is oriented through the equator at the Greenwich 
meridian, the Y E axis lies 90° to the east through the equator and the Z E axis is oriented 
up through the North Pole. The navigation function of die DAC converts these positions 
and velocities into the local Runway Coordinate System (RCS) reference frame to 
provide guidance and for post-flight evaluation. In the RCS reference frame, the origin is 
located at the aim point to the runway being flown to, the X R axis is oriented parallel to 
and down the runway, the Y R axis lies 90° to the right and the Z R axis is oriented down, 
normal to the runway. Refer to Appendix A for a full description of the RCS reference 
frame at Crows Landing NAS. 

The first step in the conversion is to transform the WGS-84 geodetic coordinates of 
the aim point being flown to into the ECEF coordinate system. Since the earth rotates, it 
assumes the shape of a sphere that is flattened at the poles and bulging at the equator. 
Therefore, the earth can be modeled by an ellipsoid of revolution formed by rotating an 
ellipse around its minor axis [10]. Figure 7 shows the ellipsoid of revolution and the 
associated ECEF reference frame. As discussed above, the ECEF coordinates (X E ,Y E ,Z E ) 
originate at the earth’s center of mass; the X E axis is oriented through the equator at the 
Greenwich meridian, the Y E axis lies 90° to the east through the equator and the Ze axis is 
oriented up through the North Pole, which coincides with the ellipsoid semi-minor axis. 
The ellipsoidal normal through a point P intersects the Ze axis, but does not pass through 
the center of mass, due to the flattening of the ellipsoid. The length of the ellipsoidal 
normal from the surface of the ellipsoid to the point P is called the geodetic height, h. 
The angle between the ellipsoidal normal and the equatorial plane is the geodetic latitude, 
<p. Any intersection of the ellipsoid of revolution with a plane containing the Z E axis is an 
ellipse called the ellipsoidal meridian. The geodetic longitude, X, is the angle between 
two meridional planes oriented clockwise from the X E axis. Therefore, the geodetic 
coordinates (<p,X,h) completely describe the position of a point on the earth. 


Zc NORTH 



Figure 7 

Earth Ellipsoid Of Revolution 


The relationship between the geodetic coordinates (<p,X,h) of the Crows Landing NAS 
Runway 35 Aim Point (AP) and the ECEF coordinates (X E ,Y E ,Z E ) is as follows [10]: 
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'ap X e ' 


(N + h)cos<pcosX, 

ap Ye 

= 

(N + h)cos<psinX, 

_ap Ze 


(N(l - e 2 ) + hjsintp 


where: 

£ 

N = , the radius of the earth ellipsoid of revolution 

^l-e 2 sin 2 cp 

e = — — - - a , the eccentricity of the earth ellipsoid of revolution 

V a a 

a is the semi-major axis of the earth ellipsoid of revolution (6378137.0 m) 
b is the semi-minor axis of the earth ellipsoid of revolution (6356752.3141 m) 
q> is the geodetic latitude of the Runway 35 Aim Point (37.41335361° N) 

A. is the geodetic longitude of the Runway 35 Aim Point (121.1082725° W) 
h is the geodetic height of the Runway 35 Aim Point (12.4 m) 

All of the numerical values above are valid only for the WGS-84 earth ellipsoid of 
revolution. 

Once the coordinates of the aim point in the ECEF reference frame are calculated, the 
difference between the airborne GPS receiver antenna position and the aim point can be 
determined as follows: 


Ax/ 


Ax e 


APxe 

Ave 

= 

Av e 

- 

ap Ye 
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- Az E _ 


APzb. 


(3) 


where: 


l x E 


L Ye 


is the position of the airborne GPS receiver antenna in the ECEF reference 


L a z e j 

frame 
AP Xe 

AP Ye 


is the position of the Runway 35 Aim Point in the ECEF reference frame 


All airborne positions, and therefore all position differences with respect to the Runway 
35 Aim Point, are based upon the position of the airborne GPS receiver antenna, which is 
the origin of the navigation solution. 

After the ECEF position difference is calculated, the airborne GPS receiver antenna 
position in the ECEF reference frame is transformed into the RCS reference frame as 
follows: 


'Axe' 


'Axe' 

Ave 

= [C"] 

Ave 

-Az R _ 


-Aze. 


( 4 ) 


where: 
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Ce = [Cv][Ce], the transformation matrix from the ECEF reference frame to the RCS 
reference frame 

-sintpcosX, -sincpsinX costp 

— sin X, cosX. 0 , the transformation matrix from the ECEF 

-costpcosX -coscpsinA, -sincp 

reference frame to the Vehicle-Carried Vertical (VCV) reference frame, where the 
origin is located at the Runway 35 Aim Point, the X v axis is oriented towards 
True North, the Y v axis is oriented towards True East and the Z v axis is oriented 
down, normal to the runway 

9 is the geodetic latitude of the Runway 35 Aim Point (37.41335361° N) 

X is the geodetic longitude of the Runway 35 Aim Point (121.1082725° W) 
h is the geodetic height of the Runway 35 Aim Point (12.4 m) 

cosH sinH 0 
-sinH cosH 0 

0 0 1 
to the RCS reference frame 

H is the Runway 35 True Heading (10.099°) 

The GPS receiver velocity is transformed, by the DAC, from the ECEF reference frame 
to the RCS reference frame in the same manner as the GPS receiver position, both 
calculated at a rate of 2 Hz, the rate at which the GPS receiver calculates its position and 
velocity. 

With the GPS receiver position solution in the RCS reference frame, approach 
guidance commands are calculated, based on the aircraft’s angular deviation from the 
nominal localizer and glideslope selected, and sent to the aircraft's approach guidance 
instruments as simulated ILS localizer and glideslope guidance at a rate of 2 Hz. 

The localizer error is calculated as follows: 

LOC Emr = tan 1 -^- ( 5 ) 

Ax* 

If the aircraft is to the right of the runway centerline, a positive localizer error is 
generated and a fly-left command is sent to the aircraft's standard approach guidance 
instruments proportional to the amount off course. Correspondingly, if the aircraft is to 
the left of the runway centerline, a negative localizer error is generated and a fly-right 
command is sent to the aircraft's standard approach guidance instruments proportional to 
the amount off course. Figure 8 illustrates the localizer error geometry. 

The glideslope error is calculated as follows: 


, the transformation matrix from the VCV reference frame 




GS Error = tan" 


f( A x,) 2 +( Ay ,) 2 


-GS s 


If the aircraft is above the selected glideslope, a positive glideslope error is generated and 
a fly-down command is sent to die aircraft's standard approach guidance instruments 
proportional to the amount off glideslope. Correspondingly, if the aircraft is below the 
selected glideslope, a negative glideslope error is generated and a fly -up command is sent 
to the aircraft's standard approach guidance instruments proportional to the amount off 
glideslope. Figure 9 illustrates the glideslope error geometry. 
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Figure 8 

Localizer Error Geometry 
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Previous flight tests involving helicopters flying steep approaches to landing 
displayed the need to provide a lateral localizer course width (full-scale lateral deflection, 
±2 dots on the Course Deviation Indicator (CDI), ) of ±350 ft at the DH [11]. In addition, 
the flight tests also indicated the need to provide vertical angular glideslope widths of 
±1°, ±2° and ±3° for the 3°, 6° and 9° glideslopes, respectively. This corresponds to a 
full-scale vertical deflection, ±2 dots on the CDI. These values were incorporated into 
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the DGPS-based guidance solution in order to give the non-standard approaches (6° and 
9° glideslopes) the sensitivity of the standard approach (3° glideslope). 

To achieve a ±350 ft lateral localizer course width at the DH, an offset aim point 
(OAP) for the localizer angular course width is utilized as the reference point for the 
angular deviations from the nominal localizer selected (See Figure 10). Since selecting a 
different glideslope changes the horizontal range (HR) between the DH and the aim point, 
a pre-defined OAP for each of the three glideslopes is utilized, depending upon which 
glideslope is selected for the approach. 
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Figure 10 

Localizer Geometry 


To achieve the vertical angular glideslope widths of ±1°, ±2° and ±3° for the 3°, 6° 
and 9° glideslopes, respectively, the aim point is utilized as the reference point for the 
angular deviations from the nominal glideslope selected (See Figure 11). Since selecting 
a different glideslope changes the vertical angular glideslope width, a pre-defined 
glideslope sensitivity (°/dot) for each of the three glideslopes is utilized, depending upon 
which glideslope is selected for the approach. 
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Figure 1 1 

Glideslope Geometry 


For all guidance and navigation calculations, a flat earth model was assumed, 
therefore, neither Coriolis or centripetal accelerations were included. In addition, gravity 
was assumed to be of a constant value. These assumptions are valid since the airspeeds 
involved are relatively slow (80 Knots Indicated Air Speed (KIAS)), and the flight tests 
are conducted at one location at relatively low altitudes (3000 ft Mean Sea Level (MSL) 
and below). 
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FLIGHT TEST PROCEDURES 


Standard (3°) and steep (6° and 9°) glideslope straight-in approaches were used to 
evaluate helicopter approaches to landing. The 3° glideslope approach was flown at 1 10 
KIAS while the 6° glideslope approach was flown at 85 KIAS and the 9° glideslope 
approach was flown at 65 KIAS. These airspeeds were chosen to keep the aircraft’s rate 
of descent less than 1000 ft/min. Fourteen separate approaches for each of the three 
glideslopes were flown for a total of 42 approaches. Throughout all of the approaches, 
the aircraft was tracked via the laser tracker. All flights were scheduled so that a 
minimum of five satellites were in view at all times and that the PDOP was less than a 
value of six. 

The overall flight path is a rectangular pattern consisting of four basic components; 
the crosswind leg, the downwind leg, the base leg and the final approach (See Figure 12). 
The climbing turn to crosswind is initiated upon reaching the departure end of the active 
runway and is flown until reaching the downwind leg. The downwind leg is offset 1 nm 
to the right of the active runway and is flown at die same altitude as the glideslope 
intercept altitude. The turn to base is initiated at 5.1 nm, which is abeam the Initial 
Approach Fix (IAF), and is the point where the data collection begins. The base leg is 
flown so that the aircraft is positioned at the IAF ready to begin the final approach. It is 
the final approach segment that is of primary interest to this research. 

All approaches were initiated at the IAF, located 5 nm out from the aim point along 
the active runway heading, with the aircraft established on speed, on course, and at the 
hard altitude associated with the glideslope to be intercepted (1 100 ft MSL, 2100 ft MSL 
and 3000 ft MSL for the 3°, 6° and 9° glideslopes respectively). Upon crossing the IAF, 
the aircraft is flown inbound to the Final Approach Fix (FAF), located 3 nm out from the 
aim point along the active runway heading. It is at the FAF that the aircraft intercepts the 
appropriate glideslope and flys the approach down to the DH corresponding to the 
glideslope just flown (190 ft MSL, 240 ft MSL and 290 ft MSL for the 3°, 6*and 9° 
glideslopes respectively). After descending through the DH, a go-around is initiated at 
which time the aircraft is flown back to the IAF to set up for another approach. A 
graphical depiction of the final approach segment is illustrated in Figure 13. 
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Figure 12 

Runway 35 Flight Test Path 
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Figure 13 

Runway 35 Final Approach 
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FLIGHT TEST RESULTS 


Data Base Summary 

Approaches to be evaluated were chosen based upon the following criteria: 1) 
continuous tracking of a minimum of four common satellites by both the airborne and 
ground-based receivers, 2) good satellite geometry (PDOP less than six) and 3) valid laser 
tracking data. 

A time history statistical analysis of the position error during the approach was the 
primary measure of the DGPS performance. Three-axis position errors for each approach 
were calculated by taking the difference between the DGPS position (in the RCS 
reference frame) with the laser tracking data (also in the RCS reference frame). In order 
to calculate the position errors, the laser tracker truth data position at the laser reflector 
was transformed to the GPS receiver antenna location as follows: 


AL Xr 

AL Yr 

=r 

Lxr 

[Cv][Cj][cJ][Cb] 

A*.' 

Av„ 

AL Zr _ 


L z ,_ 




where: 


A x .' 


'12.9540m' 

Ay. 

= 

1.4224m 

.^Za _ 


_ 3.2385m 


L the location of the right laser reflector with respect to the GPS 


receiver antenna (See Appendix A), in the aircraft body (AB) reference frame, 
where the origin is located at the aircraft’s center of gravity (CG), the X B axis is 
oriented forward along the roll axis of inertia, the Y B axis is oriented to the right 
along the pitch axis of inertia and the Z B axis is oriented down along the yaw axis 
of inertia. 


, the transformation matrix about the roll axis of inertia 


, the transformation matrix about the pitch axis of inertia 


, the transformation matrix about the yaw axis of inertia 


, the transformation matrix from the VCV reference frame 


^ W A VAVAVUVV Ui LIU 

H is the Runway 35 True Heading (10.099°) 

The GPS receiver antenna was used as the origin of the navigation solution and therefore 
all position errors are based off this location. 
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Flight tests were conducted on 11 February 1993, 2 April 1993 and 12 April 1993. 
SA was on during the flight tests. Seven 3°, eight 6° and six 9° glideslope approaches 
(out of a total of 42), which were flown on the three days with various satellite 
combinations, were chosen for the statistical analysis. The approaches that were not 
chosen for the statistical analysis had either airborne or laser tracker data that was 
incorrectly time-tagged or had incomplete laser tracker data. 

Laser Tracker Validation 

To verify the post-processing algorithm which transforms the laser tracker truth data 
to the RCS reference frame, processed data was compared to surveyed position data. 
Static calibration data was collected prior to each flight by placing the aircraft laser 
reflector over a known test point on the aircraft parking ramp (See Appendix A). Data 
collected at the test point location was transformed into the RCS reference frame in real- 
time and compared to the surveyed position. Any error greater than 0.5 m in any axis 
would require that the laser be re-calibrated and verified prior to the flight test. Laser 
calibration is verified prior to each flight test by taking measurements to several reflectors 
permanently installed at surveyed locations around the location. 

Laser Tracker Reference Results 

Figures 14 and 16 show a sample 3° glideslope approach illustrating both DGPS and 
laser tracker lateral and vertical position with respect to the AP in the RCS reference 
frame. Note the lateral position angular bias in Figure 14. This is most likely due to the 
fact that the CDI was not calibrated accurately enough to the computer-derived lateral 
guidance commands. This angular bias of approximately -0.25° was apparent in all of the 
approaches evaluated. The lateral and vertical position errors for the same sample 
approach are shown in Figures 15 and 17. The lateral position error jump at 2500 m is 
due to laser tracking ambiguities between the right and left laser reflectors (See Figure 
15). 

Composite lateral and vertical position errors are shown for all of the 3°, 6° and 9° 
glideslope approaches in Figures 18 through 23. Note that the lateral position error is 
consistently smaller than the vertical position error, as would be expected. The lateral 
position error bounds show a slight trend towards increased accuracy as the range 
decreases. Since the DGPS position accuracy is not a function of the aircraft position 
during the approach, a possible explanation would be that the laser tracker is more 
accurate at the shorter ranges. Figures 19, 21 and 23 display large vertical position error 
bounds at the beginning of the approaches. Again this is most likely due to the 
inaccuracy of the laser tracker at extended ranges. 

Note the set of data points diverging away from the zero error line at approximately 
2000 m in Figure 19. This divergence is due to the interruption of the pseudorange and 
pseudorange-rate difference correction uplink to the aircraft on flight number 3092-308 
and displays the effects of SA. These diverging data points correspond to an interruption 
that lasted for 44 sec and were not used in any position error statistical analysis. Since 
the error due to SA is a function of time, the pseudorange and pseudorange-rate 
difference corrections are time dependent. Therefore the age of the pseudorange and 
pseudorange-rate difference corrections is critical to the accuracy of the DGPS-based 
position solution. The data collected indicated that pseudorange and pseudorange-rate 
difference corrections having an age greater than approximately 20 sec resulted in 
significantly larger errors. All of the data collected was accomplished with pseudorange 
and pseudorange-rate difference corrections having an age greater than five sec, with a 
majority of the data collected having pseudorange and pseudorange-rate difference 
corrections less than 10 sec in age. 
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Figure 21 shows that the vertical position error bounds for the 6° glideslope 
approaches converge and diverge at approximately 6500 m and again at approximately 
3000 m. No correlation could be made between either the satellite geometry or age of the 
pseudorange and pseudorange-rate difference corrections and the increase followed by 
the decrease in positioning accuracy at these ranges. 

Figure 23 shows that the vertical position error bounds for the 9° glideslope 
approaches converge and diverge at approximately 2500 m. Again, no correlation could 
be made between either the satellite geometry or age of the pseudorange and 
pseudorange-rate difference corrections and the increase followed by the decrease in 
positioning accuracy at this range. 



Lateral Position (m) 


Figure 14 

Sample 3° Glideslope Approach Lateral Position 
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Lateral Error (m) 
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Figure 15 

Sample 3° Glideslope Approach Lateral Position Error 
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Vertical Position (m) 
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Figure 16 

Sample 3° Glideslope Approach Vertical Position 
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Range (m) 


Figure 17 

Sample 3° Glideslope Approach Vertical Position Error 
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Figure 18 

Composite 3° Glideslope Approach Lateral Position Error 
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Figure 19 

Composite 3° Glideslope Approach Vertical Position Error 
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Figure 20 

Composite 6° Glideslope Approach Lateral Position Error 
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Figure 21 

Composite 6° Glideslope Approach Vertical Position Error 
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Figure 22 

Composite 9° Glideslope Approach Lateral Position Error 
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Figure 23 

Composite 9° Glideslope Approach Vertical Position Error 


Time histories of the position errors were calculated at a rate of 2 Hz. A time history 
analysis assumes that the DGPS positioning accuracy is not a function of the aircraft 
position during the approach. Therefore, all of the data collected for each approach is 
given equal weighting in the statistical analysis. The composite position error statistics 
for each of the three types of approaches flown are summarized in Table 1. 


Table 1 

Composite Position Error Statistics For Each Type of Approach 


Error 

Orientation 


3° Approach Error (m) 


Mean 


SD (2o) 


6° Approach Error (m) 


Mean 


SD (2 q) 


9° Approach Error (m) 


Mean 


SD (2 o) 


Longitudinal 

Lateral 

Vertical 


-0.79 

+ 0.22 


-1.38 

-0.14 

-1.08 


-0.73 

+0.18 

-0.99 


-2.03 


±2.74 

±1.78 

±3.54 


±3.95 

±1.47 

±3.51 


±2.51 

±1.27 

±2.81 




















Since the longitudinal mean position error is biased in the opposite direction of 
aircraft motion for all of the approaches and is significantly larger than the lateral mean 
position error, it appears to be a result of computational time lags within the DGPS 
receiver, which in turn increases the age of the pseudorange and pseudorange-rate 
difference corrections. Such computational time lags were noticed at the ground-based 
DGPS receiver during the uplinking of the pseudorange and pseudorange-rate difference 
corrections to the aircraft. As the number of satellites that the ground-based DGPS 
receiver tracked increased, the rate at which the pseudorange and pseudorange-rate 
difference corrections were calculated and transmitted decreased. It should be noted that 
longitudinal mean position error is not as critical as lateral and vertical mean position 
errors during approaches to landing. The lateral mean position errors are relatively small, 
as would be expected. The vertical mean position errors are significantly larger than the 
lateral mean position errors and are primarily due to the uncertainty in the DGPS vertical 
solution, and possibly due in part to a bias due to the age of the pseudorange and 
pseudorange-rate difference corrections. 

The lateral and vertical position errors show an increase in accuracy as the glideslope 
increases (See Table 1). This is possibly due to the fact that the steeper the approach, the 
slower the airspeed and therefore the greater the number of data points to analyze, with 
the result being a more statistically accurate analysis. 

The next analysis performed was a decision height analysis, where only data points 
corresponding to the 200 ft DH on a standard 3° glideslope approach were used. This 
type of analysis is typically used in evaluating the navigational accuracy of approach and 
landing systems. The position errors at the 200 ft DH for six of the seven 3° glideslope 
approaches flown (excluding flight number 3092-308 due to the interruption of the 
pseudorange and pseudorange-rate difference corrections at the DH) are summarized in 
Table 2, as well as Figures 24 and 25. 


Table 2 

Position Error Statistics At The 3° Approach 200 ft Decision Height 


Error 

Orientation 

200 ft DH 

Error (m) 

Mean 

SD (2a) 

Longitudinal 

-0.61 

±2.79 

Lateral 

+0.34 

±1.49 

Vertical 

-2.25 

±1.63 


The results indicate improvement in position accuracy at the 200 ft DH as compared 
to the time history analysis results for the entire approach. Again this might be due to the 
increased accuracy of the laser at shorter ranges. However, it should be noted that the 
decision height analysis is based on a limited data set (six data points from six 3° 
approaches flown). 
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Figure 24 

3° Glideslope Approach Lateral Position Error Statistics 
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Figure 25 

3° Glideslope Approach Vertical Position Error Statistics 
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CONCLUSIONS 


Accuracy Considerations 

DGPS positioning accuracy is composed of lateral and vertical error components. 
The accuracy requirements for a precision approach can be expressed in terms of lateral 
and vertical error limits. A comparison of DGPS accuracy with respect to the precision 
approach requirements provides an indication of the feasibility of using DGPS to provide 
high accuracy, precision navigation and guidance for helicopter precision approaches to 
landing. 

Precision Approach Requirements For A 3° Glideslope 

The point at which the 3° glideslope is 200 ft above the surface is defined by ICAO to 
be the CAT 1 DH point [12]. The DH is defined by FAA Order 8260.3B (U.S. Standard 
for Terminal Instrument Procedures) to be the height, specified in feet above MSL, above 
the highest elevation in the touchdown zone at which a missed approach shall be initiated 
if the required visual reference has not been established. 

The lateral requirement at the CAT 1 DH point is ±17.1 m (2a) and the vertical 
precision approach accuracy requirement is 4.2 m (2a) [12]. These accuracy 
requirements are based on ICAO ILS standards for ground equipment and assume a 3° 
glideslope and 8000 ft distance between the localizer antenna and the runway threshold. 

Comparative Assessment 

DGPS positioning accuracy at the 200 ft DH on a standard 3° glideslope approach 
was 0.3 m (mean) ±1.5 m (2a) laterally and -2.3 m (mean) ±1.6 m (2a) vertically. These 
errors indicate that the helicopter position based on DGPS guidance satisfies the ICAO 
CAT 1 lateral and vertical navigational accuracy requirements (See Figures 24 and 25). 
Note that these results are based upon a limited set of data (six data points from six 
approaches). The relatively large vertical mean error is due primarily to the uncertainty 
in the DGPS vertical solution. In addition, the vertical mean error may be due in part to a 
bias due to the age of the pseudorange and pseudorange-rate difference corrections. 

Concluding Remarks 

The DGPS system utilized commercial receivers designed primarily for land 
surveying applications and was compatible with the ILS type localizer and glideslope 
instruments onboard the NASA UH-60 helicopter. Pilots commented that the DGPS- 
based guidance appeared smooth, accurate and easy to follow. 

The decision height analysis revealed that the DGPS system did achieve ICAO CAT 
1 navigational accuracy requirements for a standard 3° glideslope approach. In addition, 
14 of the 21 approaches analyzed using the time history analysis also displayed sufficient 
navigational accuracy to meet such requirements. This is a significant result, since all 
previous research, conducted with DGPS only, has not achieved such accuracy. One 
reason for this is the rapid technological advancement made in both hardware and 
software design currently available in DGPS receivers. 

The age of the pseudorange and pseudorange-rate difference corrections has a 
significant impact on the accuracy of the DGPS position solution. A bias due to the age 
of the pseudorange and pseudorange-rate difference corrections may contribute to the 
vertical position error, and to a lesser extent to the longitudinal and lateral position error. 
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RECOMMENDATIONS 


The age of the pseudorange and pseudorange-rate difference corrections needs to be 
monitored and a warning needs to be displayed to the pilot when this age becomes greater 
than 20 sec. The results of this research indicated that the positioning accuracy started to 
deteriorate significantly after a 20 sec, or more, interruption of the pseudorange and 
pseudorange-rate difference corrections. 

The guidance and navigation algorithms can be modified to increase the positioning 
accuracy. This is accomplished by integrating the DGPS and INS via Kalman filtering 
techniques. This integration is the subject of current research at the NASA Ames 
Research Center to improve helicopter precision guidance and navigation. 

In order to take advantage of the improved positioning accuracy provided by the 
DGPS/INS integration, which in turn may allow flying to a lower DH, a more 
sophisticated cockpit display unit (CDU) is needed to provide localizer and glideslope 
guidance as well as range to touchdown and overall DGPS status. Such a CDU would 
need to provide the above information in a more accurate and timely manner, as 
compared to the present instrumentation used, which would be required to fly to a lower 
DH. 
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APPENDIX A 


Crows Landing NAS Survey Data 

Previous research involving DGPS applications to fixed-winged aircraft required that 
two points be surveyed at Crows Landing NAS in 1990 [8]. These points are survey 
point AP, the aim point on Runway 35, and survey point LA on the East side of the ramp 
area (See Figure 26). The two survey points, as well as additional survey data derived 
from these two points were used in both the real-time as well as the post-processing 
portion of the current helicopter DGPS research. 

Survey point AP is located on the Runway 35 centerline approximately 3440 feet 
short of the departure end of the runway. The station mark is a brass disk stamped “AP”, 
set into the concrete of the runway, protruding up approximately 5 mm. Survey point LA 
is located approximately 6.9 meters West of the East edge of the concrete ramp (See 
Figure 27). The station mark is the center of the head of a nail sunk into the junction of 
two expansion joints of the concrete slabs forming the ramp. The expansion joints are 
oriented 350.2° magnetic (based on the current magnetic variation of 15.3° East), 
355.495° relative to the RCS reference frame and 5.594° True. Survey point L, from 
which survey point LA is defined, is located 2.579 meters South and 4.712 meters East of 
survey point LA (See Figure 27). The station mark is a brass disk stamped “L”, set into 
the concrete of the ramp, protruding up approximately 5 mm. The identification 
stamping is partially defaced. In addition, all three of the survey points are marked by 
fluorescent orange paint. 

When required to verify the laser tracker accuracy as well as the accuracy of the 
DGPS receiver or to align the INS at Crows Landing NAS, the pilot positions the aircraft 
so that the Pilot-in-Command (right seat) is directly over survey point LA. This is 
accomplished by having the pilot look over the control stick and off the right shoulder to 
align the aircraft with the expansion strips (See Figure 27). 
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Figure 26 

Crows Landing NAS Facility 
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Table 3 

Component Aircraft Body Coordinates 


LOCATION 


BUTTLINE (in) 

WATERLINE (in) 

Ground Static Line 

N/A 

N/A 

184.00 

LA 

229.00 

+24.00 

184.00 

Right Reflector 

251.00 

+56.00 

206.50 

INS 

298.75 

+32.00 

212.50 

CG 

359.50 

0.00 

258.75 

GPS Antenna 

761.00 

0.00 

334.00 


The aircraft components in the RCS reference frame are obtained as follows: 


AC Yr — LA Yr + [C B ] AACy, 


where: 

'AAC Xb ' 

AAC Yb is the difference, in the AB reference frame, of the aircraft components 
_AAC*_ 

with respect to survey point LA 
’ cosA sinA 0" 

Cg = -sin A cos A 0 , the transformation matrix from the AB reference frame to 

0 0 1 

the RCS reference frame 

A is the angular difference between the ramp orientation and the RCS reference frame 
(4.505°) (See Figure 27) 

'LA Xr 1 I" 2.5280m ' 

LA Yr = 220.0560m , the position of the point LA in the RCS reference frame 
_LA Zr J [ 0.7620m 
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Table 4 

RCS Coordinates 


LOCATION 

X R (m) 

Y R (m) 

Z R (m) 

AP 

0.0000 

0.0000 


Laser 

524.8016 

352.4311 

0.4846 

LA 

2.5820 

220.0560 

0.7620 

Right Reflector (@ Pt. LA) 

2.0348 

220.9102 

0.1905 

INS (@ Pt. LA) 

0.7778 

220.3977 

0.0381 

CG (@ Pt. LA) 

-0.8243 

219.7086 

-1.1366 

GPS Antenna (@ Pt. LA) 

-10.9909 

220.5097 

-3.0480 


Once the components are in the RCS reference frame, a transformation into the VCV 
reference frame is required in order to apply the latitude and longitude correction factors. 
This is accomplished as follows: 


AC Xv 


AC Xr 

AC Yv 

-[<S] 

ac Yr 

_ AC Zv 


ac Z r . 


(9) 


where: 
'AC Xr 


AC 

AC 


Yr 


is the aircraft components in the RCS reference frame 


Zr j 

cosH -sinH 0 


Cl = 


sinH cosH 0 


0 0 1 
to the VCV reference 


, the transformation matrix from the RCS reference frame 


rame 


H is the Runway 35 True Heading (10.099°) 

In the vicinity of survey point LA, the latitude correction factor is 1 11180.79872 m/deg 
and the longitude correction factor is 88528.5504 m/deg. Since the survey point AP can 
be used as the origin of the VCV reference frame and the latitude and longitude are 
precisely known, the component latitude and longitude are obtained as follows: 


Latitude = APla T - AC Xv /l 1 1 180.79872 m/deg (10) 

Longitude = AP^n -AC Yv / 88528.5504 m/deg (1 1) 


where: 

APlat is the geodetic latitude of the Runway 35 Aim Point (37.41335361° N) 
APlon is the geodetic longitude of the Runway 35 Aim Point (121.1082725° W) 
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Table 5 

Geodetic Coordinates 


LOCATION 

LATITUDE (N) 

LONGITUDE (W) 

ELEVATION (m) 

AP 

37° 24’ 48.073” 

121° 06’ 29.781” 

43.0 

Laser 

37° 25’ 02.830” 

121° 06’ 11.930” 

42.5 

Base GPS Antenna 

37° 25’ 04.144” 

121° 06’ 12.142” 

46.9 

LA 

37° 24’ 46.904” 

121° 06’ 20.953” 

42.2 

INS (@ Pt. LA) 

37° 24’ 46.846” 

121° 06’ 20.952” 

43.0 

GPS Antenna (@ Pt. LA) 

37° 24’ 46.471” 

121° 06’ 21.031” 

46.1 


The latitude and longitude are given in the WGS-84 reference system while the 
elevation is given in National Geodetic Vertical Datum of 1929 (NGVD 29) local 
elevation reference system. 
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APPENDIX B 


MATLAB Code For Survey Calculations 

XR, Angle between Xramp and Xr (nad) 

R-2*pi*4. 505/360; 

XLfir, Coordinates of survey point Lfl in the RCS reference frame 
(m) 

LRr“[2. 5280; 220 . 0560; . 7620] ; 

XCrampr, Transformation matrix from the ramp reference frame to 

the RCS Preference frame 

Crampr*[cos(R) sin(R) 0;-sin(A) cos(R) 0;0 0 1]; 

XGPSb, Position of the OPS antenna with respect to survey point 

LR in Xthe aircraft body reference frame (m) 

GPSb“[-1 3 . 51 28; - . 6096; -3 . 81 ] ; 

XCGb, Position of the CG with respect to survey point LR in the 
Xalrcraft body reference frame (m) 

CGb-[-3.3147;-.6096;-1 .89865]; 

XINSb, Position of the INS with respect to survey point LR in 

the Xaircraft body reference frame (m) 

I MSb-C-1 . 771 65; .2032;-. 7239]; 

XRRb, Position of the Right Reflector with respect to survey 

point LR Xin the aircraft body reference frame (m) 

RRb-[-.5588; . 81 28; - . 571 5] ; 

XH, Angle between True North and Xr (rad) 

H-2*pi*1 0.099/360; 

XCrv, Transformation matrix from the RCS reference frame to the 
Xvehi cl e-carr ied vertical reference frame 
Crv-[cos(H) -sin(H) 0;sin(H) cos(H) 0;0 0 1]; 

XRPIat, Latitude of the Runway 35 Rim Point (deg) 

RPIat«37. 4133536; 

XRPIon, Longitude of the Runway 35 Rim Point (deg) 

RPIon-121 .1082725; 

XdLRT, Latitude correction factor (m/deg) 
dLRT-1 01 .324*. 3048*3600; 

XdLON, Longitude correction factor (m/deg) 
dL0N-80. 68*. 3048*3600; 

XGPSr, Position of the GPS antenna in the RCS reference frame 
(m) 

GPSr*LRr+Crampr*GPSb 

GPSr - 

-10.9909 

220.5097 

-3.0480 

XCGr, Position of the CG in the RCS reference frame (m) 

CGr“LRr+Crampr*CGb 

CGr - 
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-0.8213 
219.7086 
-1 .1366 


XINSr, Position of the IMS in the RCS reference frame (m) 
I NSr*LRr+Crampr* I NSb 

INSr - 


0.7778 

220.3977 

0.0381 

XRRr, Position of the Right Reflector in the RCS reference frame 
(m) 

RRr«LRr+Crampr*RRb 
RRr - 


2.0318 

220.9102 

0.1905 

XGPSv, Position of the GPS antenna in the veh i c I e-carr i ed 
vertical Preference frame (m) 

GPSv-Crv*GPSr; 


XCGv , Position of 
reference Xframe (m) 

the 

CG 

i n 

the 

yeh i c I e-carr i ed 

y er t leal 

CGv“Crv*CGr; 

X 1 NSv, Pos i t i on of 

the 

INS 

i n 

the 

yeh i c 1 e-carr i ed 

yen t i ca 1 


reference Xframe (m) 

I NSv^Crv* I HSr; 

XRRv, Position of the Right Reflector in the vehicle-carried 

vertical Preference frame (m) 

RRv*Crv*RRr; 

XLRv, Position of the survey point LR in the vehicle-carried 

vertical Preference frame (m) 

LRv“Crv*Lflr ; 

PGPSIat, Latitude of the GPS antenna (deg) 

GPSIat«RPIat-(abs(GPSv(1, 1)/dLAT)) 

GPSIat - 

37.11290819696332 

XGPSIon, Longitude of the GPS antenna (deg) 

GPSI on m RP I on-(abs(GPSv(2, 1 )/dL0N)) 

GPS I on - 

121.1058120311716 

XCGIat, Latitude of the CG (deg) 

CGIat-RPIat-(abs(CGv(1 , 1 )/dLRT)) 
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CGI at 


37 .41299978550254 

XCGIon, Longitude of the CG (deg) 
CGIon«RPIon-(abs(CGu(2, 1 )/dL0M)) 

CGI on ■ 

121 .1058308020561 

XI NSIat , Latitude of the INS (deg) 
INSIat“RPIat-(abs( I NSu( 1 , 1 ) /dLflT ) ) 

I MS I at - 

37.41301 288552292 

XINSIon, Longitude of the INS (deg) 

I NS I on*RP I on-(abs( I NSu (2, 1 )/dL0N) ) 

I NS I on - 

121.1058199655112 

XRRIat, Latitude of the Right Reflector (deg) 
RRIat-RPIat-(abs(RRu(1 ,1)/dLRT)) 

RRIat - 

37.41302320790727 

XRRIon, Longitude of the Right Reflector (deg) 
RRIon«fiPlon-(abs(RRu(2, 1)/dL0N)) 

RRIon ■ 

121 .1058117769228 
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APPENDIX C 


MATLAB Code For Data Processing 

X Hn (I— file to process the rout DGP5 flight test and laser truth 
data 

GS-i nput ( 1 Enter Glideslope: '); 

RPP-input ( 'Enter Rpproach *: '); 
load RD.dat 
[r,c]=size(RD) ; 

X H, Angle between True Morth and Xr (rad) 

H=2*pi*1 0.099/360; 

X Cur, Transformation matrix from the UCU reference frame to the 
RCS 

X reference frame 

Cur-[cos(H) sin(H) Oj-sin(H) cos(H) 0;0 0 1]; 

X RLRb, Position of the right laser reflector wrt the GPS 
antenna in the 
X body reference frame (m) 

RLRb- [12. 9540; 1 . 4224; 3 . 2385] ; 

X LLRb, Position of the left laser reflector wrt the GPS antenna 
in the 

X body reference frame (m) 

LLRb-[ 1 2. 9540; -1 . 4224; 3.2385]; 
for i -I : r 

X V, Yaw angle (rad) 

Y-RD( i j 5) ; 

X P, Pitch angle (rad) 

P-RD(i,6); 

X R, Roll ang I e (rad) 

R-RD(i,7); 

X Cy, Transformation matrix about the yaw axis 
Cy«[cos(Y) -sin(Y) 0;sin(Y) cos(Y) 0 ; 0 0 1]; 

X Cp, Transformat i on matrix about the pitch axis 
Cp-[cos(P) 0 sin(P);0 1 0;-sin(P) 0 cos(P)]; 

X Cr, Transformation matrix about the roll axis 
Cr-[1 0 0 ; 0 cos(R) -sin(R);0 sin(R) cos(R)]; 

X L, Raw laser truth data 
L-[RD(i J 8);R0(i,9);RD(i JO)]; 

X RLr, Laser truth data transformed to the GPS antenna 
posit ion based 

X upon the right laser reflector (m) 

RLr-L-(Cur*Cy*Cp*Cr*RLRb) ; 

X RLI , Laser truth data transformed to the GPS antenna 
posit ion based 

X upon the left laser reflector (m) 

RL I -L-(Cur*Cy*Cp*Cr*LLRb) ; 

X Routine to determine which laser reflector is being tracked 
i f(abs(RD( i,3)-RLr(2, 1))<abs(RD(i ,3) -RLI (2, 1))) 

RL-RLr; 
e I se 
RL-fiLI; 
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end 

X D, Processed OGPS flight test and laser truth data 
D( i , 1 )-RD( i , 1 ); X IRIG-B time (msec past the top of the 

hour) 

D( i ,2)*RD( i ,2); X GPS antenna X position (m) 

D( i ,3)*RD( i ,3); X GPS antenna V position (m) 

D(i ,4) S RD( i,1); X GPS antenna Z position (m) 

D( i ,5)*RL( 1 , 1 ) ; X Laser truth X position (m) 

D( i , 6)-RL(2, 1 ) ; X Laser truth V position (m) 

D( i ,7)-RL(3, 1 ) ; X Laser truth Z position (m) 

D( i ,8)*D( i ,2)-D( i ,5); X X error (m) 

D( i ,9)-D( i ,3)-D(i ,6); X V error (m) 

D(i,10)*D(i # 4)-D(i,7);X Z error (m) 
end 

X I1EAN, Dean value of the X,Y and Z errors 
l1ERN*mean(D( : , 8 : 1 0) ) 

X TI40SIG, Two sigma standard deviation of the X,Y and Z errors 
TU0SIG*std(D( : ,8: I0))*2.0 

fprintf( 'app_errors' , 'GLIDESLOPE - X1.0f\n',GS) 

fpr i nt f ( ' app_errors ' , 1 RPPRORCH * - X-2.0f\n' ,RPP) 

fpr int f ( ' app_errors ' , ' NnXerror-mean - X-1 0 . 5f\n' , mean(D( : , 8) ) ) 

fpr i nt f (' app_errors Yerror-mean ■ X-10.5f\n‘ ,mean(D(: ,9))) 

fpr i nt f ( ' app_errors ' , ' Zerror-mean ■ X-10.5f\n\n‘ ,mean(D( : , 10))) 

fprintf( ’app-errors' , 'Xerror-2sig ■ +/-X- 

10.5f\n' , std(D( : ,8))*2.0) 

fpr int f ( ' app_errors 1 , ' Yerror-2si g ■ +/-X- 

10.5f\n' , std(D( : ,9))*2.0) 

fpr i nt f ( ' app_errors ' , 1 Zerror-2s I g ■ +/-X- 

10.5f\n\n\std(0(:,10))*2.0) 

fprintf('PD_' , 'GLIOESLOPE - XI.OfXn'.GS) 

fpr int f( ' PD_' ,' APPROACH * - X-2.0f\n' ,RPP) 

fpr int f( ' PD_' , ' \nXerror-mean ■ X-10.5f\n' , mean(D( : ,8))) 

fprintf( ’PD_' , 'Yerror-mean ■ X-1 0 . 5f\n* ,mean(D( : ,9) )) 

fprint f( 'PD_' , 'Zerror-mean ■ X-l 0.5f\n\n' , mean(D( : , 10)) ) 

fprint f( 'PD_‘ , ’Xerror-2sl g * +/-X-1 0. 5f\n ' ,std(0( : , 8))*2. 0) 

fpr int f ( ’ P0_‘ , ' Yerror-2sl g - +/-X-1 0. 5f\n ' ,std(D( : , 9))*2. 0) 

fpr I nt f( ' PD_‘ , ’ Zerror-2si g ■ +/-X-1 0. 5f\n\n' , std(D( : , 10) )*2. 0) 

fprint f( ' PD_‘ , ' IRIG-B') 

fpr I nt f( ' PD_' , ' Xgps ' ) 

fprintf('P0_', ' Ygps') 

fprintf( 'PD_' , ' Zgps') 

fprintf('PD_',' Xttl') 

fprint f( ' PD_' , ' Ytt I ' ) 

fprintf( 'PD_' , ' Zt 1 1 \n ' ) 

for 1*1 :r 

fpr int f( *PD_’ ,'X7.0f',D(i,l)) 
fpr i nt f ( ' PD_' , 'X13.5f',D(i,2)) 
fprintf( 'PD_' ,'X11.5f',D(i,3)) 
fprintf( 'PD_‘ ,'X11.5f',D(i,4)) 
fpr i nt f ( ' PD_' ,'X13.5f',D(i,5)) 
fpr i nt f ( ' PD_' , 'XI 1 . 5f ' , D( i , 6) ) 
fpr i nt f ( 'PD_' , * X1 1 . 5f\n 1 , D( i , 7)) 
end 

Xfpr i nt f ( ' PD error . dat 1 , ' IRIG-B') 
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Xfpr i nt f ( 1 PD error . dot 1 , ' 

Xfpr int f( 1 PD error . dot ’ , 1 

Xfprintf( ’PD error.dat ’ , ' 

Xfprintf( 'PD error.dat ' , ' 

for i *1 : r 

fpr intf( 'PD error.dat ' , 

fpr intf( 'PD error.dat ' , 

fpr int f ( ' PD error.dat ' , 

f pr i nt f ( ' PD error . dat ' , 

fpr i nt f ( ' PD error . dat ' , 

end 

clear; 


RANGE’) 

Xerror ' ) 
Verror ' ) 
Zerror\n' ) 


X7. Of', 0(1,1)) 
X13.5f ' ,D( i ,5)) 

X15. 1 Of ' ,D( i ,8)) 

XI 5. 1 Of ' , D( i , 9) ) 

XI 5. 1 0f\n' , D( i , 10) ) 
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APPENDIX D 


MATLAB Code For Error Analysis 

X Rn 11-file to process the DGPS Error Files 
GS“input (' Enter Glideslope: '); 

RPP-input ( 'Enter Total Number Of Approaches: 
load error.dat 
[r,c]*size(error) ; 

X flERNj dean value of the X,Y and Z errors 
l1ERN“mean( error ( : , 3 : 5) ) 

X TI40SIG, T»o sigma standard deviation of the X,Y and Z errors 
TUOSI G M std( error ( : , 3: 5) )*2 

fpr i nt f ( ' error_ana I ys i s 1 , ' GL I DESLOPE - X1.0f\n',GS) 

fpr int f ( ' error_anal ysi s 1 , ' RPPRORCHES - X-2.0f\n' ,RPP) 

fpr int f ( ' error_anal ysi s ' , ' \nXerror-mean M X- 

1 0.5f\n ' , mean ( error ( : , 3) ) ) 

fprintff 'error-analysis' , 'Yerror-mean *= X- 

10.5f\n' ,mean(error( : ,4))) 

fpr i nt f ( ' error_ana I ysi s' , ' Zerror-mean - X- 

10.5f\n\n' , mean (err or ( : ,5))) 

fpr i nt f C ’ error_ana I ys i s ' , ' Xerror-2s i g * +/- 

XI 0 . 5 f\n ' , st d( error ( : , 3) )*2) 

fpr i nt f ( ' error_ana I ys i s ' , ' Yerror-2s i g * +/-X- 
1 0 . 5f\n ' , st d( error ( : , 1 ) )*2) 

fpr int f ( ' error— ana I ysi s' j ' Zerror-2s i g ■ +/-X- 

10.5f\n\n' , std(error( : ,5))*2) 

clear; 
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APPENDIX E 


Individual Approach Time History Error Statistics 

FLIGHT NUMBER: 3042-301 

Xerror-mean ~ 0.29408 
Yerror-mean - 0.11424 
Zerror-mean - -1.90693 

Xerror-2sig - +/-1. 67565 
Yerror-2sig - +/-1. 51635 
Zerror-2sig - +/-1. 36613 

FLIGHT NUMBER: 3042-302 

Xerror-mean - -0.71309 
Yerror-mean ■» 0.29140 
Zerror-mean - -4.75735 

Xerror-2sig - +/-2. 76933 
Yerror-2sig - +/-1. 24985 
Zerror-2sig - +/-3. 42374 

FLIGHT NUMBER: 3092-301 

Xerror-mean - -2.24500 
Yerror-mean * 0.31793 
Zerror-mean • -1.36353 

Xerror-2sig • +/-1. 65130 
Yerror-2sig - +/-1. 61700 
Zerror-2sig - +/-1. 12435 

FLIGHT NUMBER: 3092-302 

Xerror-mean - -1.88567 
Yerror-mean - 0.17116 
Zerror-mean * -1.36367 

Xerror-2sig - +/-2. 49394 
Yerror-2sig - +/-1. 36208 
Zerror-2sig - +/-2. 95104 

FLIGHT NUMBER: 3092-307 

Xerror-mean - -0.23494 
Yerror-mean * 0.19723 
Zerror-mean - -0.96155 

Xerror-2sig - +/-1. 56515 
Yerror-2sig - +/-1. 53718 
Zerror-2sig - +/-0. 95922 
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FLIGHT NUMBER: 3092-308 


Xerror-mean - 0.32709 
Yerror-mean - 0.18861 
Zerror-mean - -0.57040 


Xerror-2sig - +/-1. 81597 
Yerror-2sig - +/-1. 25964 
Zerror-2sig - +/-1. 68482 


FLIGHT NUMBER: 3092-314 


Xerror-mean 

Yerror-mean 

Zerror-mean 

Xerror-2slg 

Yerror-2sig 

Zerror-2sig 


- -1.09040 

- 0.25571 

- -2.80339 

- +/-1. 71446 

- +/-3. 24581 

- +/-2. 24822 


FLIGHT NUMBER: 3042-601 

Xerror-mean — -0.50098 
Yerror-mean - -0.04846 
Zerror-mean - -1.66626 

Xerror-2sig - +/-2. 38453 
Yerror-2sig - +/-1. 50241 
Zerror-2sig - +/-1. 19721 

FLIGHT NUMBER: 3042-602 

Xerror-mean « -1.25729 
Yerror-mean - -0.02115 
Zerror-mean - -1.91473 

Xerror-2sig - +/-1. 95879 
Yerror-2sig - +/-0. 74184 
Zerror-2sig - +/-1. 48747 

FLIGHT NUMBER: 3092-603 

Xerror-mean - -3.06731 
Yerror-mean * -0.37484 
Zerror-mean « 0.16325 

Xerror-2sig - +/-1. 81117 
Yerror-2sig - +/-1. 49750 
Zerror-2sig - +/-2. 53216 

FLIGHT NUMBER: 3092-604 

Xerror-mean - -4.98386 
Yerror-mean - -0.47458 
Zerror-mean * 0.23296 

Xerror-2sig * +/-1. 83358 
Yerror-2sig * +/-1. 05925 
Zerror-2sig - +/-1. 34987 

FLIGHT NUMBER: 3092-609 

Xerror-mean * -0.90852 
Yerror-mean - 0.19866 
Zerror-mean - -1.49769 

Xerror-2sig - +/-2. 03222 
Yerror-2sig - +/-1. 86047 
Zerror-2sig - +/-6. 05498 

FLIGHT NUMBER: 3092-610 

Xerror-mean “ -0.77933 
Yerror-mean - 0.18996 
Zerror-mean - -0.56043 

Xerror-2sig - +/-1. 15870 
Yerror-2sig - +/-1. 38347 
Zerror-2sig - +/-1. 74298 
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FLIGHT NUMBER: 3092-615 


Xerror-mean - 0.19538 
Yerror-mean - -0.44835 
Zerror-mean - -2.69098 


Xerror-2sig - +/-1. 42971 
Yerror-2sig - +/-1. 58679 
Zerror-2sig ■ +/-2. 49127 


FLIGHT NUMBER: 3092-616 


Xerror-mean 

Yerror-mean 

Zerror-mean 


- 0.56622 

- - 0.20120 
0.72249 


Xerror-2sig 

Yerror-2sig 

Zerror-2sig 


+/-2. 38544 
+/-0. 85824 
+/-2. 60039 


FLIGHT NUMBER: 3042-901 


Xerror-mean - -0.69460 
Yerror-mean - -0.01253 
Zerror-mean - -1.68157 

Xerror-2sig - +/-1. 97014 
Yerror-2sig - +/-1. 03072 
Zerror-2sig * + /-1. 45544 

FLIGHT NUMBER: 3042-902 

Xerror-mean - -0.93279 
Yerror-mean * 0.17897 
Zerror-mean - -1.37699 

Xerror-2sig - +/-1. 52077 
Yerror-2sig - +/-1. 11278 
Zerror-2sig - +/-1. 99139 

FLIGHT NUMBER: 3092-905 

Xerror-mean “ -1.83832 
Yerror-mean - 0.36239 
Zerror-mean - -1.02861 

Xerror-2sig - +/-1. 93866 
Yerror-2sig - +/-1. 04739 
Zerror-2sig - +/-4. 00243 

FLIGHT NUMBER: 3092-906 


Xerror-mean * -1.65273 
Yerror-mean * 0.14191 
Zerror-mean * -2.02475 

Xerror-2sig - +/-1. 45787 
Yerror-2sig - +/-1. 68595 
Zerror-2sig * +/-2. 33202 

FLIGHT NUMBER: 3092-912 

Xerror-mean * 0.39056 
Yerror-mean - 0.24847 
Zerror-mean - -0.18022 

Xerror-2sig - +/-1. 60088 
Yerror-2sig - +/-1. 15394 
Zerror-2sig - +/-1. 55909 

FLIGHT NUMBER: 3092-917 


Xerror-mean 

Yerror-mean 

Zerror-mean 


0.22296 

0.12506 

0.12790 


Xerror-2sig 

Yerror-2sig 

Zerror-2sig 


- +/-2. 23201 

- +/-1. 35992 

- +/-1. 51378 
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